Exact low temperature results for transport properties of the interacting resonant 

level model 
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Using conformal field theory and integrability ideas, we give a full characterization of the low 
temperature regime of the anisotropic interacting resonant level (IRLM) model. We determine the 
low temperature corrections to the linear conductance exactly up to the 6th order. We show that 
the structure displays 'Coulomb deblocking' at resonance, i.e., a strong impurity- wire capacitive 
coupling enhances the conductance at low temperature. 



PACS numbers: 73.63.-b, 72.10.Fk 

Transport properties in quantum impurity problems 
have become central to experimental nanophysics [![. 
Their theoretical study involves strong interactions, often 
out of equilibrium, and thus presents considerable diffi- 
culties. In many cases, it is possible to map these systems 
to simple one dimensional models, which can be solved 
by the Bethe ansatz (BA) in equilibrium. Extending the 
BA out of equilibrium is the next step forward. This has 
been accomplished in the problem of edge state tunnel- 
ing in the fractional quantum Hall effect (where coupling 
to the reservoirs was particularly simple) B, Q , but ma- 
jor obstacles remain in most other cases [J]. Recently, 
an exciting new "open Bethe ansatz" has been proposed 
H , which might well be a breakthrough, although subtle 
issues remain open - related to universality and the treat- 
ment of boundary conditions around the impurity out of 
equilibrium: see [6] for recent progress in this direction. 

It is one of the major difficulties of this field that so few 
methods are available to investigate strongly correlated 
regimes that it is often not possible to assert the valid- 
ity of results such as those in 0] or Q (see eg 0). Our 
goal in this paper is to report on a method to tackle low 
temperature properties in the IRLM which gives highly 
non perturbative results, and apart from its own useful- 
ness, provides results that can be used as benchmarks. 
Another, different perturbative approach has been pro- 
posed recently in 

The IRLM model describes a resonant level coupled 
via tunneling junctions to two baths of spinless elec- 
trons, with which there is also a Coulomb interaction 
@> After the usual expansion in angular modes, un- 
folding and linearizing near the Fermi surface, one ends 
up with a hamiltonian H = Hq + Hb where Hq = 
— i X)a=i 2 / d% *Pl,dxip a is the free hamiltonian describ- 
ing two infinite right moving Fermi wires, and tunneling 
occurs through the impurity term: 



H B = (7i^i(0) +72^(0)) d f +h.c. 
U (V>K(0) + 4^2(0)) {J* -\)+e d 



d)d 



(1) 



In the following, it will be convenient to use the langage of 
the Kondo model, which is unitarily related to the IRLM 
@| , and to introduce spin 1/2 operators to represent the 



impurity: d' = r]S + , cftd = S z 



(77 is a Majorana 



fermion). The parameters 712 (which can be taken real) 
are parametrized as 71 + ij2 = 7v / 2e j6 '/ 2 ; 6 encodes 
anisotropy in the tunneling process. Note the presence 
of the important interaction term, a capacitive coupling 
U between the impurity and the wires. When the on- 
site chemical potential vanishes, the impurity is at 
resonance (in an actual experiment this would require 
adjusting the local grid potential to some value V g (U)). 

We shall be concerned with the conductance of the 
structure when a voltage V, that couples as Hy = 
So=i 2 I dx (ip\ip 1 — ip\ip 2 ), is applied across the impu- 
rity. The standard approach to this model is to form the 
combinations ?A + /_ = 77^75 \t 1/2^1 ^ ^2/1^2) 1 which lead 
to decoupling into two independent sectors, where H can 
be diagonalized using a straightforward BA [l(|. Using 
Friedel's sum rule the linear conductance G 
given at temperature T = by: 



dl 

dV lv=0 



IS 



G m = \ sin 



The term sin 2 6 = 2 1 1 ?\ is the familiar tunneling 

anisotropy prefactor ; the impurity filling = id' d\ can 
be extracted in 'closed form' at zero temperature using 
BA in the ip± basis However, going beyond formula 
@, which is valid only at V = T = is a very difficult 
task within this BA. The problem is that one is typically 
interested in situations where the two wires are at dif- 
ferent chemical potentials, an ensemble very difficult to 
represent in the tp± basis. We will show in this paper 
that in the field theory limit, several important results 
can nevertheless be obtained by exploiting hidden sym- 
metries. We will present the boundary conditions charac- 
terizing the zero temperature or infrared (IR) fixed point 
by making use of an underlying SU(2) structure. This 
gives an alternative, straightforward way to obtain the 
IR conductance ((2]) but also allows for the setting up of 
a convenient formalism (perturbed boundary conformal 
field theory) for systematically obtaining the low T cor- 
rections to conductance (formulae (|14)15|> ). 

Our results rely on two equivalent representations of 
([T|) in the field theory limit. We will initially set the 



sin 2 (7m d ) 



(2) 



2 



on-site chemical potential e ( i = 0, but reinstate it later. 
First, let us bosonize ib + = -7= e lv/ ^ v ±. This yields 

- 1 - V 27T 



# = Ea=± ^0(^a) + ^ [C^+ (0) S+ + h.C.] 

+ ^ (cU>+(0) + c^_(0)) 



(3) 



where the free hamiltonian is Ho(ip a ) = ^ J dx (d x (p a ) 
and k + = 7777+ . We then perform the unitary transforma- 
'(<p++v-)(o) and c hoose a = to cancel 



tion W 



Define now 



the remaining interaction along S" 2 
^[(V47r-a)( ( 9 + -Q!^_] and 0_ = ± [(y/47r- a) tp-+atp + ], 
with /3 2 = - (U — 7r) 2 + 27r. In terms of these new vari- 
ables, the hamiltonian then becomes: 

H 1 = £ Hofa) + ^± [e^+ (0) 5+ + h.c] (4) 



A remarkable feature is that the angle 9 has disappeared 
from this expression; as a result, the free energy of the im- 
purity problem is independent of the anisotropy. Hamil- 
tonian H 1 is formally equivalent to the anisotropic Kondo 
problem. The interaction term has scaling dimension 
D — f- and thus we see immediately that in the scaling 
theory, tunneling is irrelevant in the low energy limit for 
j3 2 > 8n i.e. strong Coulomb interaction \U — tt\ > ^/3tt 
and relevant otherwise. In the latter case, the 0+ hamil- 
tonian flows to the ordinary Kondo fixed point. We will 
soon argue this corresponds to the resonant level being 
hybridized with the wires, with a fixed, anisotropy de- 
pendent amount of tunneling between the two. 

To proceed, observe that we could first bosonize "01,2 
and then only form linear combinations, this time of the 
bosons. Setting ip a 
binations 0i 



^= (if! + lf 2 ), 02 = ^75 

Hamiltonian is then rotated by the same unitary trans- 
formation U = e ia ^ s * <W°) to yield 



and forming the com- 
(</?i - V2), the 



H 1 



7 



[Vi(0)O 2 (0)S + +h.c] (5) 



where we have introduced the vertex operators V±\ = 
e ±^i, V± 2 = e ±*^*», Q 2 = 7lKl V 2 + 72^2 V_ 2 and 
K a = T]T] a . The parameter (i\ = ^/2tt — a\[2 satisfies 
/3 2 = (3\ + 2ir which ensures that the perturbations have 
the same scaling dimension D. This representation of 
the Hamiltonian is more suited to non equilibrium sit- 
uation since the electrical current from wire 1 to 2 is 
simply expressible in terms of 02 only. However, H 1 
has a much more complex form than H 1 , and typically 
mixes Kondo and boundary sine-Gordon (BSG) type in- 
teractions. This is particular clear in the case j3i = 0: 
in the limit 9^0, the Hamiltonian reduces to Hg = 

2^ e iV2^4>2 S + + h c a Kondo hamiltonian. If, for the 
v 77 

sake of argument, we neglect the Klein factors (this is al- 
lowed for example in the computation of the free energy) , 
the case 71 = 72 reduces to Hg — A'j/V^tt cos(v / 27T02)5' y 
i.e. two copies of the BSG model for S y — ±\. We will 



see later that indeed, reintroducing the Klein factors, the 
model at /3i = and 71 = 72 shares the same (Dirichlet) 
boundary conditions with the BSG model. 

That the free energy of the two incarnations H 1 and 
H 1 is the same and independent of 9 is remarkable. It 
can be checked directly but not straightforwardly at all 
order of the perturbative expansion in powers of 7. 

We consider the question of the linear conductance at 
low energy, i.e. when the resonant level is hybridized with 
the wires. A quick way to obtain it is to use the boundary 
conditions (BC's) for the fields 0i,2 in the IR. These are 
not so easy to obtain from hamiltonian H 1 . We will thus 
start from the BC's for the fields 0±, which are known 
from the general analysis of the Kondo model. The idea 
is to follow these BC's through the canonical transfor- 
mations. This seems very hard due to the non linearities 
involved, but becomes possible once one recognizes the 
presence of SU(2) affine currents (the SU(2) transforma- 
tions are those mixing the two wires). Introducing 



J* = ^sin(x/8^0 2 ) , jy = ^ C0S (V8^>2), 
(here p = k 1 k 2 ) it is easy to show that 



(6) 



d x (j)- 



7T~ dx<i 



2tt^ [cos (9 J 2 +sin6LP] 



_ 



f d^i + ?f [cos 9J Z + sin 61 J*] 







(7) 



We now recall that in the IR the field 0+ obeys Neumann 
BC's with angle f , + (O+) = 0+(O~) + §, and the boson 
0_ - being unaffected by the interaction - Neumann BC's 
with angle 0. Introducing SU(2) rotated currents J° = 
TZg-J a (TZg is a rotation of angle 9 around J v ), we see that 
in the IR, these currents obey the BC's J z (0+) = J Z (CT) 
and ^(0+) = -J±(Cr). A little algebra based on the 
SU(2) commutation relations then leads to 



J z (0+) = cos(260 J Z (Q-) +sm(29)J x (CT) 
J x {0 + ) = -cos(20) J*(0~) + sin(26>) J 2 (CT 
J v (0+) = -J y (Q-) 



(8) 



which of course are highly non linear in terms of the field 
4>2 itself. As for 0i, it obeys simply 0i(O + ) = 0i(O _ )+^ L . 
The BC's for 02 interpolate continuously between Neu- 
mann (N) (6» = and 2 (O+) = 2 ((r) + v /V8) and dou- 
ble Dirichlet (D) (0=§ and 2 (O+) = -0 2 ((T)± 
This is possible because the dimension of the operator 
e ±i\/27r02 j s the inverse of an integer square, here ^ flij ]. 
As a result, the ratio of boundary degeneracies for N and 
D is g N /g = 2, and degeneracies of N and double D 
are the same. That the two fixed points can be reached 
depending on 9 is particular clear in the case f3\ = dis- 
cussed previously, where H M interpolates between Kondo 
and double BSG. Independently of 9, the IR boundary 
degeneracy thus always takes the N value, g 1R = g N . 

Another way of viewing these BC's is to observe that 
the radius of compactification r = v2 being an integer 



3 



multiple of the self dual radius r* — allows for the ex- 
istence of a pair of non chiral operators of dimension one, 
e ±i v / 27r(^2+02)^ w hj c h induce an exactly marginal bound- 
ary deformation through J x -y- z . To our knowledge, this 
is the first time such a BC is encountered in a condensed 
matter context. It is also the first time that a tunable 
flow into a line of boundary fixed points is encountered. 

A similar analysis can be carried out when there is a 
chemical potential for the electron on the dot, i.e. a term 
eddld = €d(S z + \)- As before we first argue for the 
Kondo hamiltonian H 1 . The local magnetic field results 
in an additional phase shift for the electrons [13j, which 



translates into a phase shift 5+ for the field 



This is 



easy to understand: an additional phase shift for </> + in 
the IR is tuned by the introduction of a scattering po- 
tential term, which then induces, by Friedel sum rule, 
an extra "charge" (here, magnetization) on the impurity. 
Now, it is known that the problem with a field acting on 
the impurity only is closely related [14( (in the scaling 
limit) to the problem with a field coupling to the total 
spin S z ot = S z + jj J d x (/)+ which is immediately solv- 
able by BA. This gives rise to the impurity magnetization 



nij 



(S z 



n d - 



i, with the relation 5 4 



At T = 0, the impurity magnetization can be obtained 
using the Wiener Hopf technique. One finds two possible 
expansions. For small enough one has 



1 



(-l) r 



W 2ti+1 ^ 
1 \2(1-D), 



''imp 



,271+1 



DJU ^ n\(2n + 1) W (2 " +1)Z3 1 



(9) 



with 



u = -7= T r(j/2(i-D)) rt ' ^ e ^ ua ^ ex P ress i° n 
should be used beyond the radius of convergence u* = 



1 



in, 



20? 



E 



(-l)T(l+n(l-J)) ^ 
nl T(l - nD) 



(10) 



The parameter Tb is the Kondo temperature for the 
problem defined by hamiltonian H 1 , related to the bare 

coupling by Tb/W = r -^§^E§f^ (^ , 
with W the bandwidth (this relation holds within the reg- 
ularization inherited from integrability used in Ref.[lj|). 

The additional phase shift for 4>+ translates into more 
complicated non linear BC's for the boson <p2 '■ namely, 
the rotated currents J a now obey the BC's J z (0 + ) = 
J Z (CT) and ^(0+) = -e ±lA J ± (0-), with the angle A = 
^nTftifrip . .After some calculations one finds 

J z (0+ ) = [cos 2 9 - sin 2 9 cos A] J z (CP ) 

+ sin(26>) cos 2 f J x ((r) + sin sin A J v (CP ) 

J*(0+) = [sin 2 9 - cos 2 9 cos A] J x (0~) 

+ sm(20) cos 2 f J z {0-) -cos(9sinA Jf((T) ^ 

J y (0+) = - sin 9 sin A J 2 (CT) 

+ cos 9 sin A J x (0" ) - cos A JV (0~ ) 

These BC's relate the currents on both sides of the im- 
purity through a SU(2) rotation now depending on the 
anisotropy and doping of the impurity. 



To extract information from the BC's in the IR, it is 
convenient to reformulate them first within a boundary 
field theory by folding and introducing complex coordi- 
nates z = t — ix, x > 0. The Kubo formula then reads: 



0= I'm -i- / dxdydre 1 ^ (j e (x,T)j e (y,0)) (12) 



c^o h (2L) 2 u 



where < t < T _1 , spacial integrals run over [0, L] and 
L has to be sent to +00; the electrical current through 
the whole structure is j e {x) = 2( y J z (x) — J z (~x)) (the 
conductance depends only on the 4>2 propagator). 

Using the propagators that can be deduced from (| 1 lj) : 



4tt {d x 4>2{z)d x (f> 2 {w)) = (z-w) 
4tt {d x (t>2{z)d x 4> 2 {w*)) = (l - 2 sin 2 6>cos : 



(13) 



one finds G IR = %- sin 2 9 cos 2 — , which is nothing but 
Eq.@. The capacitive coupling U has disappeared at 
the IR fixed point: it is 'irrelevant', but as we will see, it 
controls the approach of the fixed point and thus deter- 
mines the low T < Tb properties of the theory, which we 
will now be able to tackle, thanks to this long reformu- 
lation of Friedel sum rule. 

Indeed, the exact solution of the Kondo hamiltonian 
leads to a full knowledge of the infinity of counter-terms 
necessary to describe the approach to this fixed point 
[l5j ]. allowing one to carry out IR perturbation theory 
to all orders. A program such as the one of [3| could 
then lead to results for the linear conductance at arbi- 
trary values of the temperature. It relies on the identi- 
fication of the low T Hamiltonian, which has the form: 
H = H IR + J2k>0^2k-iO2k{x — 0). It is important to 
stress that this expansion is highly non perturbative in 
the tunneling amplitude (as we will see below, it leads to 
an expansion of the conductance in powers of 7~ 2 / 1_£l ). 
But the couplings &2fc-i turn out to be known explic- 
itly [l(| ! The whole set of perturbing operators C>2k is a 
set of commuting conserved quantities related to integra- 
bility and describes the approach to the IR fixed point; 
it is made of fields of even dimensions. 2 k can be ex- 
pressed as a polynomial in d x <fi+ and its derivatives, to be 
then translated in the 01,2 basis. We just sketch here the 
(somewhat lenghty) analysis. Apart from density-density 
couplings, the leading irrelevant contribution contains a 
tunneling term 0| un — -M'0l'0i+'02' ! /'2)(' ! /'iV'2+ n - c -) w ^ n 
a coupling constant A cx Ipsinfl. The anisotropy and 
coulombic repulsion are only apparent in the amplitude 
of the tunneling term - which, as it should, vanishes in 
the Kondo limit, 9 = 0. This pattern generalizes to all or- 
ders: the whole set of operators describing the approach 
to the IR fixed point is independent of U. 

To obtain the conductance, the current-current corre- 
lator in Eq. (TT2")) is expanded in powers of the couplings 
&2fc-i- The resulting multiple integrals over intermedi- 
ate times of finite T correlators are evaluated using the 
residue theorem; divergences are regularized in the 'in- 
tegrable' scheme through the commutativity of the 02fc- 



4 




FIG. 1: First three reduced coefficients a 2k = g2k(TB/^) 2k 
(defined in (|14[0 . The scaling dimension D varies between j 
and 1, which corresponds to the region \U/n-l\ < V3 where 
tunelling is relevant. The inset displays the universal ratio 



9o94/ 



which diverges for D 



(U 



This way, we extract the low T expansion of G, yielding 
the 'Landau-Fermi parameters' gik for the conductance: 



G = G 1R (l + ]T.9 2fc T 2fc ). 



(14) 



fc>0 



It is important to realize that the coefficients in this ex- 
pansion are universal in the field theory limit. They can 

be put in the form g 2k — a>2k\T^) with a 2fc depending 
only on U. Laborious calculations yield: 



AX 



0,2 



04 



3(1 + Xf 
16A 
45(1 + X) 2 



(15) 



3X 



{i + xy 



X{lh-X) 

' 167T 



where we introduced the parameter X — AD — 1, and y = 

° n figE the fl2fc ' S are pl ° tted 
up to order 2k — 6 (we have obtained 06 as well, but its 

expression is too lengthy to be shown here). 

The lowest order correction g 2 can be understood in 
a simple way: at generic values of the coulombic re- 
pulsion on the impurity, the IR fixed point is a Fermi 
liquid whose approach is controlled at lowest order by 
a single operator, namely the energy momentum tensor 
(d(j) + ) 2 . Now corrections to G IR can only stem from this 
part C>2 Un of the perturbing operator that involve charge 
transfer across the impurity site, whose amplitude is sim- 
ply multiplicativcly rcnormalized with respect to the free 
case: X(U) = fqrj|A(0). This reasoning fails for higher 
orders, that are controlled by several processes with dif- 
ferent couplings, each of them being a function of U. 

Note that the £/-dependance of Tg results in a maxi- 
mum for G in the scaling limit ^/VW <g; 1 ; this effect, 
which comes from the boundary perturbation being most 
relevant for D = j (corresponding to U = 7r in our renor- 
malization scheme), was also noted in at T = 0, V ^ 
using perturbation theory in U. 

The Landau-Fermi parameters allow to form a number 
of universal ratio, the simplest one being p = = 2|. 
The first order of its development in U agrees with re- 
sults in 0. It displays (see fig.l) a divergence at the par- 
ticular value U = it of the Coulombic repulsion, which 



might offer an efficient way to identify this point in ex- 
perimental realizations of the IRLM. Moreover, at this 
value of U, there are no processes allowing for charge 
transfer up to order 6: the coefficients a 2 and 04 van- 
ish - indicating the somewhat singular nature, for trans- 
port properties, of the Fermi liquid at this point - while 

This results in a further cn- 



«6 



2tt 



' 105 Vr(2/3) a , 



hancement of the conductance around U = tt (see fig. 2) 
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FIG. 2: Iso-conductance plot in the (U/W,j 2 /W) plane, at 



fixed T/W=0.0005. On each line, G/Gm = 1 
ranging from 1 (dark) to 7 (bright). 



10~ x with x 



In conclusion, it should be clear that methods of field 
theory give one a complete control of the linear conduc- 
tance problem from the IR point of view. Apart from 
their practical use (the 8 th order could be calculated and 
the series Pade resummed to obtain full crossover curves) 
we hope that our results will provide useful benchmarks 
in assessing other approaches to the problem. 

The authors would like to thank B. Doyon and A. Kom- 
nik for numerous enlightening discussions, as well as N. 
Andrei and P. Mehta for their patient explanations. 
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